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\Q ; ABSTRACT 

■ We present a catalogue of positions, magnitudes and velocities for 3300 emission-line objects 



found by the Planetary Nebula Spectrograph in a survey of the Andromeda Galaxy, M3 1 . Of 
Q^, these objects, 2615 are found likely to be planetary nebulae (PNe) associated with M31. The 

I ' survey area covers the whole of M3 1 's disk out to a radius of 1 .5°. Beyond this radius, obser- 

O \ vations have been made along the major and minor axes, and the Northern Spur and Southern 

Stream regions. The calibrated data have been checked for internal consistency and compared 
^ ' with other catalogues. With the exception of the very central, high surface brightness region 

' of M31, this survey is complete to a magnitude limit of msoo? ^ 23.75, 3.5 magnitudes into 

the planetary nebula luminosity function. 

We have identified emission-line objects associated with M3rs satellites and other back- 
rS , ground galaxies. We have examined the data from the region tentatively identified as a new 

■ satellite galaxy, Andromeda VIII, comparing it to data in the other quadrants of the galaxy. We 

find that the PNe in this region have velocities that appear to be consistent with membership 
ofM31 itself. 

The luminosity function of the surveyed PNe is well matched to the usual smooth mono- 
tonic function. The only significant spatial variation in the luminosity function occurs in the 
vicinity of M3 1 's molecular ring, where the luminosities of PNe on the near side of the galaxy 
are systematically ^0.2 magnitudes fainter than those on the far side. This difference can be 
explained naturally by a modest amount of obscuration by the ring. The absence of any dif- 
ference in luminosity function between bulge and disk suggests that the sample of PNe is 
not strongly populated by objects whose progenitors are more massive stars. This conclusion 
is reinforced by the excellent agreement between the number counts of PNe and the R-band 
light. 

The number counts of kinematically-selected PNe also allow us to probe the stellar dis- 
tribution in M31 down to very faint limits. There is no indication of a cut-off in M3rs disk 
out to beyond four scale-lengths, and no signs of a spheroidal halo population in excess of the 
bulge out to 10 effective bulge radii. 

We have also carried out a preliminary analysis of the kinematics of the surveyed PNe. 
The mean streaming velocity of the M3 1 disk PNe is found to show a significant asymmetric 
drift out to large radii. Their velocity dispersion, although initially declining with radius, flat- 
tens out to a constant value in the outer parts of the galaxy. There are no indications that the 
disk velocity dispersion varies with PN luminosity, once again implying that the progenitors 
of PNe of all magnitudes form a relatively homogeneous old population. The dispersion pro- 
© 2006 RAS ^jjj asymmetric drift results are shown to be mutually consistent, but require that the disk 

flares with radius if the shape of its velocity ellipsoid remains invariant. 
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1 INTRODUCTION 

Observations of the distribution of stellar velocities in the solar 
neighbourhood of the Milky Way reveal a wealth of complex struc- 
ture (.Pehnen & Binnev...l998.) . For a dynamical entity such as the 
Galaxy, this structure is as fundamental as the spatial arrangement 
of its stars, and should contain a mass of information as to the sys- 
tem's current state, as well as archaeological clues as to how it 
formed. Unfortunately, such data from a single locality are difficult 
to interpret in terms of the global structure of the galaxy. Although 
kinematic observations of mo re distant stars in the M ilky Way have 
been made (see, for example. ISevenster et aljfl997h . the complex- 
ity of disentangling the geometry of our own galaxy make such data 
difficult to interpret. We also clearly need more than one example 
before any general conclusions about the dynamics of disk galaxies 
can be drawn. An obvious next target is therefore the Milky Way's 
nearby sister system, the Andromeda Galaxy, M31. 

However, studying the detailed kinematics of stars in even a 
nearby external galaxy is quite challenging. In fact, the first prob- 
lem arises from its very nearness: M31 subtends such a large angle 
on the sky that many spectroscopic instruments are unable to survey 
the whole system within a realistic time. Further, even at what by 
extragalactic standards is th e very modest distance of 785 ± 25 kpc 
iMcCoimachie et al]l2005l) . obtaining high-quality spectra of in- 
dividual stars is challenging; it is only relatively recently that the 
requisite observations of e ven bright i^d giant branchJRGB) stars 
have been made feeitzel & Guhathakurtall2002t llbata et alJl20o4 
iKalirai et &\VM) ^. 

Planetary nebulae (PNe) have long been recognised as a po- 
tentially simpler dynamical tracer of the stellar population. Stars 
evolve very rapidly but quite gently from the RGB phase to become 
PNe, so the kinematics of these two populations should be essen- 
tially identical. Further, the presence of strong emission lines in 
PNe make them quite easy to identify, and render the measurement 
of their Doppler shifts fa irly trivial. More than twenty years ago, 
iNolthenius & Ford h984) reported on kinematic measurements of 
34 PNe in M31, but it was not until recently with the development 
of wide-field multi-object spectrographs that it has become pos- 
sible to start to obtain th e large sample necessary for dynamical 
studies jHurlev-Keller et a l. 2004; Hallidav et al. 2006). Even with 
such instruments, these studies involve both narrow-band imaging 
to identify candidate PNe and subsequent spectral follow-up, and 
this complex two-stage process has limited astronomers' enthusi- 
asm for such projects. However, the recently-commissioned Plane- 
tary Nebula Spectrograph has reduced this process to a single stage, 
simultaneously identifying PNe and measuring their velocities. We 
have therefore used this instrument to carry out a deep kinematic 
survey of PNe over ~ 6s square degrees, covering the bulk of 
M3rs visible area. 

This paper describes the survey, tabulates the final data set of 
3300 detected emission-line objects (of which 2730 are probably 
PNe, 2615 in M31), and presents some initial analysis of these data. 
In Section |2| we briefly describe the instrument and the observa- 
tions made. A detailed description of the data analysis is presented 
in Section |3| Section |4| compares the photometric and kinematic 
results with those from previous smaller surveys, while Section |5| 
explores the contamination of the data by objects other than PNe 
and the sample's completeness. The catalogue resulting from all 
this analysis is presented in Section |5| In Section we identify 
PNe associated with M3 1 's satellites and other galaxies in the field, 
both to study these systems and to create a "clean" M31 sample. We 
end the paper with a preliminary analysis of the resulting sample of 



M31 PNe: Section [8| investigates their luminosity function; Sec- 
tion |5| compares their spatial distribution with that of the starlight 
in the galaxy; Section fTol calculates the luminosity specific number 
density of PNe; and Section fTTI examines the basic kinematic prop- 
erties of M3 1 's disk out to large radii that can be gleaned from these 
data. We draw conclusions and look to the future in Section fT2l 



2 INSTRUMENTATION AND OBSERVATIONS 

2.1 The Planetary Nebula Spectrograph 

The Planetary Nebul a Spectrograph (PN.S), described in detail in 
[Dou 2las et al. ( 2002), was designed specifically to measure the 
shifted wavelength (and hence velocity) of the 5007 A [O III] emis- 
sion line that dominates the spectra of PNe. It makes this mea- 
surement using a technique called "counter-dispersed imaging." 
Light from the galaxy initially passes through a narrow-band fil- 
ter 35 A wide) centred on the emission line. The light is then 
divided into two beams, with each beam dispersed by a grating be- 
fore being imaged onto separate CCDs, termed the left- and right- 
arm images. As illustrated in FigureQ we thus obtain a pair of slit- 
less spectrographic images in which stellar continua appear as elon- 
gated streaks, or "star trails", and PNe are visible as dots from their 
[Qui] emission lines. However, the two gratings are arranged to 
disperse the light in opposite directions, so the dots are offset from 
the true location of the PN in opposite directions as well. Thus, in 
essence, the mean location of the dots in the two images gives the 
position of a PN, while their separation yields its velocity. 

The instrument's ability to obtain all the necessary infor- 
mation from a single observation, combined with its relatively 
large field of view,^ make it a very efficient tool for study- 
ing PN kinematics in external galaxies. Its main drawback rel- 
ative to the conventional approach of optical identification and 
follow-up spectroscopy is that the attainable velocity accuracy 
is somewhat lower at ~ 15-20 kms~^, compared with th e ~2- 
10 k ms~^ obtain ed from traditional spectroscopy ( Hallidav et ajJ 
.2006; iHurle^Keller et al. 2004; CiarduUo et al. 2004) . However, 
this level of accuracy is good enough to study even quite subtle 
kinematic properties of external galaxies. 

2.2 Observations 

The PN.S observations for this survey were carried out on the 
William Herschel Telescope (WHT) in La Palma on two runs: 2002 
October 8-13 and 2003 September 29 - October 5. In M31's halo, 
where PNe are scarce, supplementary imaging was used to optimise 
the choice of fields to be observed with PN.S. This imaging was ob- 
tained using the Wide Field Camera on the Isaac Newton Telescope 
(INT), also in La Palma, on 2003 August 16 - 21. Figure|2|shows 
the locations of the fields observed with both instruments. 

Fifteen minute exposures were used for the majority of PN.S 
fields, although multiple observations were made for a few fields. 
INT WFC exposures were twenty minutes long through the narrow- 
band [O III] filter and two minutes with the off-band g' filter. For 
the PN.S observations, we used the instrument's 'A' filter at a 0° 

^ Mounted on the William Herschel Telescope using the Isaac Newton 
Group's EEV CCDs, the PN.S has a field of 10f39 by 11.'35, with pixel 
scales of 0'.'301 per pixel in the dispersed direction and 0'.'272 per pixel 
in the non-dispersed direction, and a spectral dispersion of 0.775 A per 
13.5 |am pixel. 
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Figure 1. A sample of PN.S data. The left hand panel shows a conventional [O III] image of a small section of M31 iMassev et all200^ . while the two right- 
hand panels show the same field as seen through the two arms of PN.S. Note how continuum objects are dispersed into "star trails" while [O III] emission-line 
objects remain as shifted point sources. 
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Figure 2. Fields observed in M31. The small square fields are PN.S field locations, with numbers 1-163 observed on the first PN.S ran, the rest on the second; 
the larger fields with dotted outlines are the INT WFC fields. The ellipse marks a 2° (27.4 kpc) disk radius. 
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Figure 3. Seeing as measured for different survey fields. Darker squares 
represent better seeing, with values in arcseconds as given on the accompa- 
nying scale. Empty squares contained no detected PNe. 



3 DATA ANALYSIS 

3.1 Wide Field Camera data reduction 

The INT Wide Field Camera data frames were debiased and 
flat-fielded using standard data reduction packages in IRAF. The 
g'-band images were remapped to match the coordina tes of the 
[O III] images, and SExtractor jSertin & Arnout3ll99d) was used 
in double-image mode to find sources in the [O III] image and si- 
multaneously measure the flux of counterparts in the g band. Ob- 
jects with more than five times as many [O III] counts as g' counts 
were selected as potential PN candidates for PN.S follow-up. These 
candidates were then visually confirmed individually, and their ap- 
proximate positions were calculated; more accurate astrometry was 
not required as the PN.S has a large field of view. 

3.2 Planetary Nebula Spectrograph data reduction 

Since the PN.S is such an unconventional instrument, we have had 
to produce a customised data reduction package to convert the raw 
data to a scientifically-usable form. This package has been con- 
structed within IRAF, supplementing the standard tasks with For- 
tran and PGPLOT extensions. Given the non-standard nature of this 
procedure, we now describe the reduction package and its applica- 
tion to the M3 1 observations in some detail. 



tilt angle JPouglas et all2002ft . This configuration provides a cen- 
tral wavelength of 5002.2 A and a 36.5 A bandwidth, allowing 
PNe to be detected with velocities between —1369 kms^^and 
817 kms-\ This range is well matched to the velocity range 
of components of M31, which lie between ~ —600 kms~^and 
km s^^, but also allows a significant margin for the detection of 
any unexpected high- velocity features. 

The adopted survey strategy was to tile the entire area of 
M3rs disk out to a 2 degree radius with overlapped PN.S pointings. 
We also fully mapped regions of particular interest at larger radii, 
such as the Northern Spur and Southern Stream (Ibata et al. 2001). 
However, at larger radii PNe become sufficiently scarce that most 
fields would be completely empty, so we used the INT narrow-band 
imaging to identify candidate PNe, and only made PN.S observa- 
tions of fields where at least one object had been detected. In total, 
226 PN.S fields were observed. 

Weather conditions varied considerably during the PN.S runs 
with two half nights lost to cloud cover during the first run, and 
four nights lost in the second. The seeing also changed significantly 
over the runs, with values between O'.'S and 3'.' 1 recorded. To make 
optimum use of these varying conditions, we observed the central 
fields in the worst seeing conditions (see Figure |3j: PNe are suf- 
ficiently plentiful in these fields that we did not need to get so far 
down the luminosity function to obtain a useful dynamical sample, 
and the high continuum surface-brightness in this region meant that 
we would in any case not be able to achieve completeness at fainter 
magnitudes (see Section|5}. Confusion with discrete objects is not 
a significant problem in surveys of this type as stars must be fairly 
bright to leave a visible trail. Even the very brightest stars at the dis- 
tance of M3 1 only leave faint stellar trails and in the bright central 
regions they are not visible above the background light. 



3.2.1 Initial pipeline 

Images were first debiased by subtracting a surface function fit to 
the pre- and over-scan regions using the 7i?AFtask imsurSt. There 
was very little structure in the bias of the EEV CCDs, so the func- 
tion subtracted was close to flat. 

To map out and remove bad regions on the CCDs, zero- 
exposure bias frames were used to create a map of pixels where 
charge transfer problems have occurred. Skyflats were employed 
to locate other low-sensitivity pixels (by comparing to a median 
smoothed version; pixels with greater than 9a difference between 
the two were taken to be bad). These maps were then combined 
with a list of known bad pixel regions to make a mask which is 
fed to the standard IRAF routine fixpix to correct all the flats and 
science images. 

Normalised pixel response maps were made by dividing flats 
by a median-filtered image (to remove spatial structures) then com- 
bining the flats using weights (to keep the noise Poissonian) and 
rejection (to eliminate cosmic ray events). The science and calibra- 
tion images were flatfielded by dividing through by this response 
map. 

The next step is to remove any residual cosmic rays: if they 
were left in through the next stages that remap pixels, then they 
would become smeared out, and could well mimic the point-spread 
function of a PN, so it is important to eliminate them at this stage. 
The cosmic rays were identified and remove d using t he Laplacian 
cosmic ray identification routine, lacosjm ivan Dok kum 200 ij), 
with five iterations and a contrast limit between the cosmic rays 
and the underlying object of 1.2. This produces fairly clean images 
and any residual cosmic ray events are unlikely to be a problem 
as there would have to be a pair of residuals in the two images in 
locations such that they could be a PN. 

The resulting cleaned science images still contain curvatures 
and rotations which must be eliminated in order to match up the 
left-right image pairs. The distortions involved are somewhat com- 
plicated for a slitless spectrograph like the PN.S, since they involve 
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Figure 4. Detailed arc spectrum seen through filter-A set at 0° . The wave- 
lengths of all the detected lines within the bandpass are annotated. 



both imaging distortions and the wavelength calibration of the dis- 
persed light. To determine the requisite mapping, we obtained fre- 
quent calibration images in which a mask containing 178 regularly- 
sp aced holes is introd uced into the beam (as described in detail 
m iDouglas et ai]l2002h . By illuminating this mask with a copper- 
neon-argon arc lamp, we obtained an array of short filter-limited 
spectra on the CCD, allowing both the spatial mapping and the 
wavelength calibration of the instrument to be determined. As Fig- 
ure |4| illustrates, even the small wavelength range admitted by the 
narrow-band filter contains a number of lines. We used the brightest 
of these lines (at 5017 A) to measure any departures in the locations 
of these dots from the regular array of holes in the mask, thus deter- 
mining the spatial mapping. These distortions were then removed 
from the accompanying science images using IRAF's standard geo- 
metrical mapping routines, geomap and geotran. Two other lines at 
4990 A and 5009 A were then used in conjunction with the 5017 A 
line to characterise the wavelength calibration across the image by 
interpolating a quadratic solution between these points. A database 
of these solutions for the array of points on each CCD was pro- 
duced; by interpolating between adjacent solutions, we can map 
directly from the locations of an object detected in the two arms, 
{xL,yL,XR,yR} to a true position on the sky and an emission 
wavelength, {xo,yo, A}. 

Following the spatial correction step, it is also possible to co- 
add images in those cases where multiple observations of the same 
field have been taken. Any offsets between images were determined 
by measuring the positions of dispers ed brighter stars in the field 
(selected from the USNO-B catalogue ^Monet et all2003h . and the 
necessary shifts were made. As a slight refinement, we actually per- 
formed both the image undistortion and this shift in a single trans- 
formation of the original image, so that there is only one interpo- 
lation of the data. We then combined the matched image frames 
using a weight determined by the quality of each image. These 
weights were calculated by using the IRAF routine Gtprofs to fit 
a one-dimensional Gaussian across each star trail used in determin- 
ing the shifts; the parameters of these fits were used to quantify the 
transparency, seeing and background level of each exposure. The 
appropriate weighting for each exposure was then given by 
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Figure 5. Difference in undispersed y coordinates of matched pairs of M3 1 
point sources detected in the two arms of the PN.S. The dotted line shows 
those few pairs deemed spurious by the application of further cuts at Yr — 
Yl = -1.80 and 1.17. 



3.2.2 Identification of emission-line objects 

Emission-line objects in the pipeline-processed images were iden- 
tified by a semi-automated routine. The undistorted images from 
both left and right arm were median-subtracted to remove unre- 
solved continuum from M31, and SExtractor was run on the image 
to detect point sources. In order to eliminate the incorrect identi- 
fication of features in star trails as point sources, SExtractor was 
rerun on the same image after it had been convolved with an ellip- 
tical Gaussian function elongated in the dispersion direction; any 
sources originating from the stellar spectra were still found as this 
blurring has essentially no effect on star trails, so these false detec- 
tions can be eliminated. 

Sources independently identified from the two arms of the 
spectrograph were initially automatically matched in left-right 
pairs by selecting detections that have the same coordinate in the 
undispersed y directions to within two pixels, and whose coordi- 
nates in the dispersed x direction differ by an amount consistent 
with the bandpass of the filter and the wavelength calibration de- 
termined above. These candidate pairs were passed on to a custom- 
written Fortran routine which allows the two images to be examined 
side-by-side. Any false detections, such as confused points within 
extended structures, were eliminated at this point. Further manual 
inspection then added back in any detections that have been missed 
by the automated routine, usually due to proximity to the field edge 
or a star. The edited list of emission-line pairs was then passed to 
the JRAF phot routine to improve the position determination and to 
measure the brightness of the source [using an aperture radius equal 
to the mean full width at half maximum (FWHM) seeing for that 
field to minimise sky contamination]. The FWHM of each source 
was determined using the Moffat parameter returned by the IRAF 
imexam task, as this turned out to be the simplest robust metric. 

Once the source positions on both arms have been refined in 
this way, we can apply a slightly more stringent cut in matching up 
the spatial coordinates of sources to eliminate chance alignments: 
as Figure|5|shows, the distribution of differences in the y coordinate 
between the two arms is significantly tighter than the initial cut of 
±2 pixels, so we deemed the objects in tails of the distribution to 
be chance alignments, and eliminated them from the source list. 

At this point, we have the definitive list of pairs of objects, 
and by applying the calibration determined in Section l3.2. II we 
can transform these coordinates into a spatial location and a wave- 
length. Identifying the emission with the 5007 A [O III] line, we 
can translate the observed wavelength into a velocity, which we 
correct to a heliocentric value on a field-by-field basis using JRAFs 
rvcorrect routine. A measure of the internal consistency of these 
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Figure 6. Differences between the velocities obtained for repeated observa- 
tions of PNe in the overlaps between fields. A Gaussian fit gives a dispersion 
fromPN.S of 17 kms-^. 



velocities can be obtained using the PNe that lie in the overlaps 
between fields for which we have multiple velocity measurements. 
There are 732 objects detected in more than one field (643 in two 
fields, 88 in three fields, 1 in four fields), so we have quite a large 
sample to play with. Figure|S|shows the distribution of velocity dif- 
ferences between pairs of measurements. A Gaussian fit to these 
data gives a dispersion of 24 kms~^, implying an error on each 
individual measurement of 17 kms"'^. This value is clearly an 
over-simplification since the distribution appears somewhat non- 
Gaussian and possibly even multi-modal. This complex shape can 
be attributed to the complexity involved in the wavelength cali- 
bration of slitless spectroscopy. However, the repeat observations 
used in Fig. |6| all lie close to the edges of the instrument's field 
of view, where distortions in the wavelength solution are greatest 
and the amount of calibration data is smallest. The quoted error 
should therefore be viewed as conservative, but with some system- 
atic residual errors on the scale of ~ 5 — 10 kms^^ not ruled out. 
Nonetheless, even these pessimistic errors are entirely adequate in 
a study of the large-scale kinematics of a massive galaxy like M31. 

3.2.3 Astrometry 

The measured coordinates must next be transformed into a stan- 
dard reference frame. This process was carried out using the same 
basic routines used for image stacking (see Section l3.2.U . The ap- 
proximate location of each field was determined from its nominal 
pointing, and star positions in the vicinity were read in from the 
USNO-B catalogue. The corresponding stars trails in the PN.S im- 
ages were found and locations calculated using IRAF?, xregister 
task; the resulting coordinates were fed to the ccmap task to calcu- 
late the astrometric solution for the field. In the very central fields, 
the bright background means there are very few stars available. 
In these fields, FN positions that have already been calibrated in 
overlapping neighbouring fields have been added to the list of co- 
ordinates, in order to bootstrap the astrometric solution across the 
central region. 

The PNe in the overlap regions between fields also provide an 
internal check on the accuracy of the astrometry, by measuring off- 
sets between the coordinates for the same source in different fields. 
As with the velocity consistency check described above, this will 
tend to be a conservative measure of accuracy, since spatial dis- 
tortions are greatest near the edges of fields where overlaps exist, 
so one might expect the astrometric solution to be poorest in these 
regions. As Figure0shows, the astrometric accuracy is somewhat 
better in declination than right ascension. This difference arises be- 
cause right ascension corresponds to the direction in which the cal- 
ibrating stars are dispersed, so centroiding their emission in this di- 
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Figure 7. Differences between the positions obtained for repeated obser- 
vations of PNe in the overlaps between fields. Gaussian fits to these give 
measurement uncertainties of o-ra = 0"34 and ctdcc = 0''16. 



rection is fundamentally less accurate. However, the absolute posi- 
tional accuracy is more than adequate for this large-scale kinematic 
survey of M3 1 . 

The astrometric matching of fields also allows us to make a 
final pass through the list of detected sources to eliminate the du- 
plicate detections from field overlaps. Where the duplicate observa- 
tions were obtained under similar seeing conditions, the positions, 
fluxes and velocities of the two detections were simply averaged; 
where the seeing conditions vary by more than (X'S, just the better- 
quality data were used. 

For this survey, it is also useful to obtain the coordinates in an 
M31 -based reference fra me. We have therefor e followed the geo- 
metric transformations of IHuchra et alJ il99lh . which define right 
ascension and declination offsets relative to the centre of M31, 
{£,, rj}, and coordinates aligned with the system's major and minor 
axes, {X, Y}, via the relations 



^ = sin(RA - RAo) cos(Dec), 

rj — sin(Dec) cos(Deco) 

— cos(RA — RAo) cos(Dec) sin(Deco) 

X = -^sin(PA) - 77Cos(PA), 

Y = -^cos(PA) +r?sin(PA). 



(2) 

(3) 
(4) 
(5) 



here, we have taken the centre of M3I to lie at RAo = 
00'^42™44^3, DecQ = 41°16'09" (J2000.0), and have adopted a 
positi on angle for the major axis of PA — 37? 7 Ids Vaucouleurd 
Il958h . With this choice of coordinates, positive X is located south- 
west of the centre of M3 1 and positive Y lies to the northwest. 



3.2.4 Flux calibration 

Although not primarily designed as a photometric instrument, the 
PN.S does provide information on the magnitudes of detected PNe 
through the brightness of the spots in the two arms. However, the 
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Table 1. Instnimental efficiencies calculated from spectrophotometric stan- 
dard stars. 



Star 




Date 




Efficiency 










Left 


Rieht 


Total 


LDS 7496° 


154.759 


2002-10-10 


0.1424 


0.1467 


0.2891 


BD-i-33 2642'' 


6367.221 


2002-10-10 


0.1305 


0.1322 


0.2628 


/"> 1 ni 1 AC 


63.438 


ZUUz-iU-i J 


0.1533 


0.1393 


0.2926 


BD+17 4708= 


18303.018 


2003-09-29 


0.1454 


0.1435 


0.2890 


BD-i-28 4211 = 


8708.624 


2002-10-08 


0.1365 


0.1325 


0.2691 


BD-l-28 4211 = 


8708.624 


2002-10-10 


0.1331 


0.1290 


0.2621 


BD-l-28 4211 = 


8708.624 


2002-10-11 


0.1328 


0.1299 


0.2626 


BD-l-28 4211 = 


8708.624 


2002-10-11 


0.1322 


0.1312 


0.2634 


BD-l-28 4211 = 


8708.624 


2002-10-12 


0.1331 


0.1298 


0.2629 


BD-l-28 4211 = 


8708.624 


2003-09-29 


0.1438 


0.1414 


0.2852 


BD-l-28 4211 = 


8708.624 


2003-09-30 


0.1371 


0.1292 


0.2662 



References: '' lo3<1974 . MOke & Gun3<1983h . = lok3ll990t) 



Table 2. Extinction per unit airmass in r' from the Carlsberg Meridian Tele- 
scope, and weather notes for the nights observed. 



Night 


Ext^/ 


EiTor 


Extsooo" 


EiTor 


Comments 


2002-10-08 


0.101 


0.004 


0.172 


0.007 


Rain at end 


2002-10-09 


0.129: 


0.029 


0.219: 


0.049 


Cloud first half 


2002-10-10 


0.106 


0.006 


0.180 


0.010 




2002-10-11 


0.107 


0.005 


0.182 


0.009 




2002-10-12 


0.106 


0.005 


0.180 


0.009 




2002-10-13 


0.112 


0.009 


0.190 


0.015 




2003-09-29 


0.234: 


0.024 


0.398: 


0.041 




2003-09-30 


0.213: 


0.012 


0.362: 


0.020 


Cloudy at start 


2003-10-01 


0.317: 


0.068 


0.539: 


0.116 


Mostly cloudy 


2003-10-02 










Mostly rain 


2003-10-03 










Rain and cloud 


2003-10-04 










Rain 


2003-10-05 










Rain 



" The conversion to Extsooo is calculated from TableQin RGO/La Palma 
Technical Note 3 1 to be a factor 1 .70. 5000 A is the closest listed value to 
the filter central wavelength. 

conversion of count rates into fluxes is not entirely straightforward. 
The PN.S has a rather slow shutter that takes some 15 seconds to 
fully open or close, so the effective vignetting in short exposures 
is significantly different from that in long exposures. This limita- 
tion means that we could not use sky flats to determine the overall 
vignetting, so instead used long exposure dome flats obtained dur- 
ing the day. To obtain an optimal vignetting function, dome flats 
obtained at many position angles were averaged to render the il- 
lumination of the field as uniform as possible. Since our science 
exposures are 15 minutes long, we have not applied any correction 
for the small variation in effective exposure time due to the slow 
shutter. 

We measured the total system efficiency for PN.S in the obser- 
vational setup described using observations of spectrophotometric 
standard stars. The flux per unit wavelength was measured from a 
narrow column through the brightest region of the stellar trail and 
compared with the published value. As Table Q shows, the total 
system efficiency estimates are found to be surprisingly consistent, 
with a mean value of Efftotai = 27.3 ± 0.4%, with much of the 
scatter attributable to the non-photometric conditions. As only a 
few standard stars observations were made during the run this level 




Am, 



(mog) 



Figure 8. Magnitude variations within the PN.S data set from PNe in 
field overlaps. A Gaussian fit to this gives a measurement uncertainty of 
^^"15007 = 0.07 mag. 



of consistency may be artificially low. However, since the science 
goals of this kinematic survey do not require high precision pho- 
tometry, this uncertainty is of little significance. 

Fluxes of the detected sources were then calibrated using the 
formula 

_ {CL.gL + CR.gR)E5 

Efftotal 



5007 



'5007 



(6) 



where Cl and Cr are the total counts from the left and right arms; 
ql and qr are the gains of the CCDs on each arm; E5007 is the en- 
ergy of a photon at 5007 A; A is the effective geometric collecting 
area of the telescope and is equal to 13.85 m^; Toxp is the exposure 
time; asky is the nightly sky extinction as recorded by the Carls- 
berg Meridian Telescope (see Table |2j multiplied by the airmass 
at which t he field was ob served; and acai the Galactic extinction 
taken from lSchleeel et alj (199^) and corrected to our filter central 
wavelength using the relation from lCardelli et all ^98^ giving a 
value of OGai ~ 0.226. The conversion to PN-specific magnitudes 
follows the Jacobv 1 1989) relationship. 



7715007 = -2.51og(F5007) " 13.74, 



(7) 



where the zero-point is chosen such that an emission line object 
would have the same apparent magnitude as would be recorded if 
observed through a V-band filter. To give a sense of this zero-point 
for objects in our catalogue, the resulting magnitudes range from 
a fairly bright H II region at 17.8 mag to the faintest PNe at ~ 26 
mag. 

Once again, we can obtain an internal estimate of the uncer- 
tainty in magnitudes by comparing the values derived for dupli- 
cate observations in the overlaps between fields. Figure |8| shows 
the dilferences between magnitudes in these repeated observations; 
a Gaussian fit to this distribution gives a rather small measurement 
uncertainty of Um^oor ~ 0.07 mag. Magnitudes for extended ob- 
jects will be systematically underestimated by this technique as no 
allowance has been made in the choice of aperture size for objects 
that are larger than the field point-spread function. However, since 
the main focus of this survey is to determine the properties of un- 
resolved PNe, this issue is also not a major concern. 



4 COMPARISON WITH OTHER DATA SETS 

One concern with any data set from a novel instrument like the 
PN.S is that it may contain unpredicted systematic errors. Fortu- 
nately, a number of smaller data sets obtained using more conven- 
tional instrumentation already exist, so we can compare the new 
data with these subsets to check for such systematic effects. 
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Figure 9. Comparison of initial PN.S velocities to iHallidav et alJ I2O 

fibre-spectrograph velocities. The solid line shows a fit to the data, whUe 
the dotted line is at kms~^. 




Figure 12. Plot of the difference between the final PN.S and the 
iHurlev-Kell er et alJ l2004) velocities. A gaussian fit gives At) = —5 ± 
1 kms^^ and ct = 12 ± 1 kms~^. 




'PNS 'Fibre 



Figure 10. P l ot of the difference between final PN.S velocities and 
IHallidav et alj i200d) fibre-spectrograph velocities. A gaussian fit gives 
A?; = -0.7 ± 0.4 kms-i and a = 15.4 ± 0.4 kms"!. 



4.1 Velocity comparison 

The largest existing kinematic survey of PNe in M3 1 is that by 
[Hallidav et al. (2006, hereafter H06), which determined velocities 
for 723 PNe using the conventional approach of narrow-band imag- 
ing and fibre-fed spectroscopy. Of these PNe, we find 715 objects 
within 4"of a PN.S object; the remaining 8 are likely obscured be- 
hind stellar trails in the PN.S images, but a 99% recovery rate is 
very respectable, and already gives a measure of the completeness 
for these brighter PNe. 

If we plot the difference between H06 and PN.S velocities 
against the PN.S velocity, as in Figure|9| we find a small zero-point 
offset and a very slight linear dependence on velocity, fitted by the 
relation 



WPN.S — UH06 



: 0.0151DPN.S - 17.0. 



(8) 



This small systematic effect presumably arises from the limi- 
tations of the restricted range of arc lines available for PN.S wave- 
length calibration (see Figure |4}. There is no large scale spatial 
dependence in this velocity difference, though there appears to be 
some small variation within the PN.S field. It makes very little dif- 
ference to any of the scientific results, but for consistency we have 
applied the correction implicit in equation Jsj to our complete data 
set and for all subsequent analysis. The residual differences in ve- 
locity between PN.S and H06 data are shown in Figure lTol Fitting 
a Gaussian to this distribution gives a dispersion of 15.4 kms^^; 
H06 quote an error of 6 km s^^ for their fibre spectroscopy, imply- 
ing that the PN.S data have an uncertainty of 14 km s~ ^ . This value 
is somewhat smaller than the error deduced from self-calibration 
for the reasons discussed above, and is probably a more accurate 
measure of the true uncertainty in the PN.S velocities. 

An estimate of M31's system velocity for the PN.S data was 



made by averaging the mean velocities from radial bins on either 
side of the major axis for data within 1° (~ 2.3 scale lengths) of the 
centre so as to exclude warps and asymmetries that have been ob- 
served at large radial distances, leading to a value of —309 kms^^. 
This measurement is in r easonable agreement with the standard 
value of — 300±4 km s~^ ide Vaucouleurs et alJl99lh and in good 
agreement with H06's value of —309 kms"'^. 

As a test of the robustness of this velocity cross-calibration, 
we can also compare our results with a sm aller but differently lo- 
cated survey carried out by IHurlev-Keller et alj 12004). They sur- 
veyed one quadrant of M3rs halo, finding a total of 146 emission- 
line objects, 135 PNe associated with M31, eight in M32 and three 
H II regions in Andromeda IV (see Figure II H . We have covered 
the majority of their survey area, finding 126 of their emission- 
line objects. The distribution of differences in velocities is shown 
in Figure fT2l A Gaussian fit to this distribution gives a mean ve- 
locity difference of — 5 ± 1 kms^^ and a combined dispersion of 
12 kms^^, lower than the claimed PN.S velocity error. However, 
the distribution of errors here is clearly non-Gaussian, so not too 
much weight should be given to the exact value. Nonetheless, the 
errors are clearly small, well below a level that would compromise 
this kinematic study of M3 1 . 



4.2 Astrometric comparison 

The PN.S astrometry can also be compared with the H06 data set. 
The distribution of differences in coordinates are shown in Fig- 
ure fT3l The right ascension distribution is similar to that found 
internally within the PN.S data, with a combined dispersion of 
~ 0'.'6. In declination there is a small systematic offset of -l-0'.'l8, 
with a long tail to positive values. As Figure fTsl shows, the tails 
of the distribution are almost entirely populated by PNe from the 
four central PN.S fields. As discussed above, the lack of measur- 
able stars in these fields compromises the PN.S astrometric solution 
somewhat, so such larger errors are not surprising. Excluding these 
fields, the combined dispersion in declination is ~ 0'.'2. Given the 
enormous size of M3 1 , such uncertainties are completely negligible 
for kinematic purposes. 



4.3 Photometric comparison 

Although not principally a photometric survey, we have been able 
to obtain fairly consistent magnitudes from the PN.S data, so once 
again it would be useful to compare against other data to obtain an 
external measure of the data quality. Fortunately, M3 1 has been ex- 
tensively imaged as part of the Local Group Survey i Massev et alJ 
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Figure 11. Positions of other surveys with respect to the PN.S sur vey. The PN.S survey area is shown in grey; Local Group Survey fields iMasse^^JlQQ^ 
are the large fields with solid outlines; the lCiardullo et al] ^1989^ data are the solid black points near the centre; the open circles are the lHurlev-Keller et alj 
|2004|) data; and the small points are H06's fibre spectroscopy data. 



I2OO2I hereafter M02), which includes narrow band studies of both 
the [Olll] line and Ha + [Nll] lines. As Figure fTTI shows. M02 
imaged 10 fields over the majority of the disk of M31 and hence 
a large portion of the PN.S survey region (covering 2745 of the 
objects in the PN.S survey). 

We therefore analysed the M02 survey images by obtaining 
fluxes at the positions of PN.S sources using the IRAF phot task 
with a fixed aperture of 2"- a value somewhat larger than the typ- 
ical survey seeing (~ O'.'S — I'.'S) to allow for positional uncer- 
tainties, but not so large as to introduce significant source confu- 
sion. The photon counts were converted to fluxes using factors of 
3.92 X 10"^'^ and 1.79 x 10"^'' for the [O III] and H a -1- [N II] data 
respectively (Massey, private communi cation). Galactic extin ction 
corre ctions were made using data from lSchleeel et all il998h and 
ICard elli et al. ( 1989), as for the PN.S data, yielding extinctions of 
0.225 mag and 0.164 mag for [Olll] and H a -1- [N II] respectively. 

Figure ^1 shows the comparison between the imaging mag- 
nitudes derived in this way and the PN.S data. If we remove the 
faintest PNe for which photometry is rather uncertain, and proba- 
ble H II regions for which this aperture photometry is inappropri- 



ate, there is a good overall agreement between the data sets. Fitting 
a Gaussian to the remaining 1769 PNe gives negligible zero-point 
offset of 0.02 mag, and a combined dispersion of 0.16 mag. There 
is no evidence for varying offsets in different fields or for data taken 
on different nights. 



As a further cross-check, Figure fTsl shows the comparison be- 
tween the PN.S magnitudes and the smaller data set published by 
Ciardullo et al. 1 1989). As expected, the scatter increases somewhat 
for the fainter PNe, but the overall agreement is very good. A Gaus- 
sian fit to the magnitude differences yields a mean of 0.00 mag with 
a combined dispersion of 0.13 mag, and even at the faintest mag- 
nitude the combined dispersion is less than 0.3 mag. Once again, 
there is no evidence of varying offsets for different fields or data 
taken on different nights. Since we are only going to use the pho- 
tometric data to make quite crude cuts in the luminosity function 
of PNe, which we have measured over a range of more than three 
magnitudes, photometry with an error of less than 0.3 magnitudes 
is all that is required. 
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Figure 13. Compaiison between PN.S and H06 astrometry, with differences 
measured in arcseconds. The black Hne represents all the PNe found in both 
catalogues; the grey shaded area shows those PNe in the central PN.S fields 
(60,61,70,71) 
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Figure 15. Comparison between PN.S magnitudes and those published in 
Ciardullo et al. 1 1989). Examining this plot by eye there appears to be a 
slight trend, with PN.S magnitudes being fainter than the corresponding 
[ciardullo et al. 1 1989) magnitudes for faint PNe. However, gaussian fits in 
half magnitude bins show this to be small (< 0.1 inag) and of low signifi- 
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Figure 14. Photometric compaiis on between the magn itudes derived from 
the PN.S data and those from the lMassev et alJ i2002h imaging data. The 
black outline represents all the objects that are within both surveys. In the 
upper panel the shaded area represents PN remaining after a number of cuts 
have been applied to the data. The lower panels show the distribution of 
objects for the cuts indicated. The larger, light grey area is the total cut 
made while the smaller dark grey histograms show the particular objects 
being excluded by each criterion. 



fw/ seeing 

Figure 16. A plot of the FWHM of sources detected by the PN.S survey, 
normalised by the seeing of each observed field. Sources in the dotted long 
tail toward large sizes have been flagged as contaminating H 11 regions. The 
cut between extended and point sources is made at fw/seeing = 1.25. 



5 CONTAMINATION AND COMPLETENESS 

The next issue to be addressed is to try to identify those emission 
hne sources that are not PNe, and to estimate what fraction of PNe 
we may have missed in the survey. 

With only a single emission hne detected by PN.S, any object 
that emits [O III] may form part of the data set (or, indeed, any ob- 
ject along the line of sight with a redshift that places an unrelated 
emission line at this wavelength). The dominant contaminating ob- 
jects are likely to be Hll regions, which emit strongly in [Olll], 
and exist in large numbers in disk galaxies like M3 1 . Most of these 
objects will already have been excluded, since their extended na- 
ture means that they will have been cut from the data set at the 
source-identification stage (even the largest PNe, with diameters of 
up to 3 pc, will be essentially unresolved in these data). However, 
H II regions have a broad size distribution, and the more compact 
ones will be barely resolved so will still be in the sample. Figure fT3l 
shows the distribution of source sizes, normalised to the seeing, of 
the objects remaining in the survey; there clearly remains a tail of 
resolved objects, so we have flagged sources for which the FWHM 
normalised by seeing exceeds a value of 1.25 as probable H II re- 
gions in the final catalogue. 

When examining the spatial distribution of the objects flagged 
by this process, it was found that they largely lie in a ring around 
MSl's centre. A further clump of 35 objects was also found in a 
single field close to the centre of M3 1 . This field was the only one 
close to the galaxy's centre that was observed under good seeing 
conditions (~0'.'8; see Figure |3}- With the seeing this good, the 
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Figure 17. The flux ratio of [OlII] to (HQ + [Nn]) versus magnitude. 
The upper panel shows all the data, with extended objects shown as black 
crosses. The lower panel shows the same data as histograms in different 
magnitude bins, with the extended sources shaded. The vertical line in each 
histogram shows a suggested division between PNe and H II regions based 
on these distiibutions. These lines have been plotted as triangles in the upper 
panel, and the dashed line connects them. 



variation in instrumental focus across the field and differences be- 
tween the two PN.S images may become significant, so it is possi- 
ble that these apparently extended objects are actually unresolved 
PNe. However, for consistency we also flag them as non-PNe; as 
discussed previously, we in any case do not expect the survey to be 
complete in these crowded high surface-brightness fields. 

Presumably, this cut on angular size will still miss the most 
compact H II regions and in the central regions where the seeing 
was very poor objects with diameters of up to ~ 15 pc will be 
classified as compact by this cut. Fortunately, we do have some 
further information that we can use to try to identify these re- 
maining contaminants. Specifically, we have measured fluxes from 
both the [Olll] and (Ha -I- [Nil]) lines from the M02 data, and 
the ratio of these fluxes, R, differs between PNe and H II regions 
iCiarduUo et al. 2004). In particular, since PNe give out most of 
their light in the [O III] line, we would expect R to be large for 
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Figure 18. Locations in position and velocity of sources identified as proba- 
ble H II regions. The whole sample is shown as small dots, with the probable 
H IT regions as larger crosses. 



these objects. Attempts have been made to quantify this criterion, 
leading to the su ggestion that f or bright PNe one should place a 
cut at i? ~ 1 - 2 iciarduUo et~l. 2002, 2004; Maerini et al. 2000). 
However, it became clear when we started investigating this quan- 
tity for the data in the PN.S survey that the appropriate value for this 
cut varies significantly with magnitude. As Fi gure fTTl illustrate s . the 
optimum cut in R to exclude H II regions (as calibrated using the 
identified extended sources) decreases at fainter magnitudes. We 
have therefore applied the cut shown as a dashed line in Figure [TTl 
with all sources below the line flagged as H II regions, regardless 
of their sizes. Although the distribution of extended sources on 
this plot would indicate that some fraction of the sources above 
thi s line are H II regions and th is c ut is lower than the o nes made 
bv lCiarduUo et alJj200l l2004') and' Magrini et aljj200 0>). we have 
been chosen to be conservative by only selecting for exclusion the 
area of the plot that seems virtually devoid of PNe. 

Figure fTsl picks out the locations in the and {X,v} 

planes of the objects that these two criteria have identified as non- 
PNe. Although the criteria were designed to flag up probable H II 
regions, it is clear that they also exclude most of the objects with 
velocities that are inconsistent with membership of M31. Closer 
inspection shows that the majority of these sources have emission 
lines extended in the dispersion direction, indicative of Doppler 
broadening of the emission from massive background galaxies. We 
could, for example, be observing the redshifted A3727 A [Oil] 
emission from star-forming galaxies at redshifts ^ 0.34. Clearly, 
more extensive conventional spectroscopy would be required to 
confirm any such identifications. Figure [Tsl also shows the struc- 
ture in the distribution of H II regions associated with spiral arms 
as well as the ring mentioned above. At a radius of 0.8°, the ring 
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Figure 19. Cumulative distribution of radii in the disk plane for bright 
(solid line), intermediate-luminosity (dotted line) and faint (dashed line) 
PNe. A Kolmogorov-Smimov test shows the distribution of bright and 
intermediate-luminosity PNe to be statistically indistinguishable at the 95% 
confidence level. There is a dearth of faint PNe near to the centre of M31 
resulting from poor seeing conditions and a high surface brightness. 



coinci des with the star-forming ring discussed by IPevereux et all 
il994) . where an excess of H II regions might be expected. 

As a final check on our reliabihty in identifying contaminating 
sources, we have compared our identifications with the ca talogue of 
1312 emission-line objects that lMevssonnier et alJil993l hereafter 
MLA93) detected and classified in M31 using broader-band slitless 
spectroscopy. We find 856 of these objects to be within 4" of an ob- 
ject in our survey, the great majority of which are classified as PN 
by both MLA93 and ourselves. Only 94 of these objects are listed 
by MLA93 as non-PNe, and we have successfully flagged two- 
thirds of these as non-PNe as well. Of the remaining 31, MLA93 
identify 25 as possible Wolf-Rayet stars. Since Wolf-Rayet stars 
are at a late stage in stellar evolution rather similar to PNe, and the 
two are sometimes even found together iGornv & Stasiiiska 1995), 
it is not clear that they should be considered as contaminants at 
all. One object classified as a possible QSO or Wolf-Rayet star is 
discussed further in Section The last five objects are of uncer- 
tain classification in the MLA93 catalogue and will be left in this 
catalogue. The 456 objects that were not found presumably do not 
exhibit [O III] emission. 

The other side to this coin is to ask how many PNe we might 
have missed. A simple test of completeness can be made by com- 
paring the distribution of PNe with radius in different magnitude 
bins. As Figure [T9l shows. bright and intermediate-luminosity PNe 
show spatial distributions that are statistically indistinguishable, but 
the faint PNe display a dearth of objects at small radii. As discussed 
above, such incompleteness is expected given a combination of the 
high surface brightness of the bulge region and the poor seeing con- 
ditions during observations of the central fields. Overall, we can 
tentatively conclude that the sample is complete to msoo? = 23 
over the entire survey, to m^oor ~ 23.75 at projected disk radii 
larger than 0.2°, and to m^oor ~ 25 beyond 1°; we will revisit this 
issue when we look at the PN luminosity function in Section|8| In 
terms of spatial coverage, as Figure |2| shows, M31 has been com- 
pletely mapped out to a deprojected disk radius of 1.5°. Beyond 
that radius, the major and minor axes are somewhat unevenly sam- 
pled from one side to the other, but still with reasonable spatial 
coverage. Thus, although not spatially complete, useful kinematic 
information can be gleaned out to 2° (27.4 kpc). 



6 THE CATALOGUE 

Having calibrated the data spectroscopically, astrometrically and 
photometrically, both internally and against external results, we are 
now in a position to present the complete database from the sur- 
vey. The full final version of the resulting data table of 3300 emis- 
sion line objects including 2730 probable PNe is included in the 
electronic version of this paper, but Table |3| gives a small sample 
of the tabulated data listing coordinates, magnitudes, heliocentric 
velocities, whether the source is a probable PN, whether it lies 
in M31 (see Section Q, when the data were obtained, and cross- 
identifications with other catalogues. 



7 SATELLITE GALAXIES AND OTHER OBJECTS IN 
THE PN.S SURVEY AREA 

Within the area of this survey are a number of galaxies besides 
M31. Some of the objects in our catalogue which reside in these 
systems can be identified in the database because their velocities are 
far from that of M3 1 . As we have seen above, most such sources 
have already been flagged as non-PNe because of their extended 
nature or low [Oiii]/Hq ratio, and are probably in quite distant 
faint background systems. In some cases, such as the sources in the 
less-distant background galaxy Andromeda IV, there may be a few 
genuine PNe. Of more immediate interest are the satellite galaxies 
around M3 1 . In these systems, we do detect quite a number of PNe, 
which means we can obtain some crude measure of the internal 
kinematics of these systems, but we must take some care in flagging 
these sources so that they do not compromise the use of this data 
set to model the dynamics of M3 1 . 

The locations of the known galaxies in the survey field are 
indicated in Figure l20l and we now consider them in turn, to try 
to identify the sources that lie within them, both on the basis of 
their positions and their velocities. Probable non-members of M3I 
identified in this way are annotated as such in the full database (see 
Section|6}. 



7.1 M32 

M32 (NGC 221), is a compact dwarf elliptical galaxy (cE2) pro- 
jected in line with the disk of M3 1 . Its major and minor axes mea- 
sure ~ 8^7 and 6.'75 in diamete r respectively to a B-band magni- 



tude level of 25.0 mag arcsec Ide Vaucouleurs et all 1 991 ). It has 



a system velocity of —200 ± 6 km s ^ iHuchra et alll99S ), and a 
velocity dispersion of 50 ± 10 km s~^, rising to ~ 80 kms~^ at 
its centre ( Simien & Prugniei2002h . 

We have considered all the emission-line objects in the PN.S 
survey within a 5^7 radius (1.5 times the mean 25.0 mag arcsec"^ 
B-band radius) of the centre of M32 as candidate members of the 
satellite (see Figure 1211 . Using the PNe at somewhat larger dis- 
tances from M32 to characterise the kinematics of M3rs disk in 
this region, we find that a velocity cut at —350 km s~^ (three times 
the measured dispersion from M32's systemic velocity) yields a 
reasonably clean division between the two systems. However, it is 
possible that a few of the PNe on the M31 side of the cut may be 
objects with the most extreme velocities in M32. To warn of this 
possibility, objects in the catalogue are flagged as "M32" if they 
meet the cut in velocity as well as position, and "M32?" if they are 
simply coincident with M32's location. 

The sample of likely M32 PNe has only 46 members, so there 
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Figure 21. PN.S data in the vicinity of: (a) M32, (b) NGC 205, (c) And IV, (d) 2MASXi J0039374+420956, (e) MLA93 0953. Tiie area around each minor 
galaxy is shown in the upper panels. The sizes of each symbol indicate an emission-Hne object's velocity with respect to M31, with crosses showing negative 
values and stars showing positive values. The large circles indicate the approximate spatial extent of M32 and NGC 205 (1.5 times the mean 25.0 B-mag 
arcsec"^ radius) and the location of the smaller galaxies. For M32 and NGC 205 PNe with circles around them have both position and velocity consistent 
with membership of that galaxy. The lower panels shows velocity histograms, the dotted line being the whole PN.S sample, the solid line indicating an area 
including the minor galaxy but extending well beyond its optical extent, the hatched histograms showing the PNe within the area of galaxy, and the cross- 
hatching highlighting PNe that agree in velocity as well as position; the points indicate the smaller galaxy's systemic velocity (from the literature, if known), 
with 3(Tu errorbars. 
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Table 3. Emission line objects identified in tlie M3 1 survey. See electronic version for full table. 



ID 


RA 


Dec 


"^5007 


^hclio 


Note 


Host 


Date 




Field 


H06 C89 


HK04 MLA93 




J2000.0 


J2000.0 






a 


h 








Number 


c c 


c c 


1 
i 


U.41.Ui. / 


4Z.0 /.jo.O 


23.48 


-2o3.3 






ZUUZ 


1 
lU 


1 9 
12 


1 






2 


U.4U.JD.4 


42. 1U.32.U 


22. lo 


-133.9 






9AA9 
ZUUZ 


1 A 
lU 


1 A 
10 


2 






i 




4z.oy.3y.D 


25.51 


O/iO A 

2oo.y 




2MAOA1 


9009 
ZUUZ 


1 
lU 


1 A 
10 


3 






4 


(J. jy.j /.3 


A O .AO. CO *i 

4z.oy.jy.o 


23.52 


3j4.j 


E 


2MA^>A1 


9AA9 
ZUUZ 


1 
iU 


1 A 
10 


3 






5 


u.jy.ji.o 


4z. 1 1 .JO.O 


23.28 


-343.3 






9AA9 
ZUUZ 


1 A 
lU 


10 


3 






6 


0:38:54.0 


42:09:39.3 


25.36 


-266.7 


E 




9009 
ZUUZ 


1 
lU 


10 


4 






1 


0. j /.j4.3 


/IT- 1 /I -/IQ Q 
4Z. 14.4y.3 


23.3 / 


no 1 
-32o. 1 






9009 
ZUUZ 


1 
lU 


1 A 

10 


c 

J 






8 


0:43: 1 1.6 


A O .A"7. 1 1 A 

42.0/.13.U 


22.31 


-194.0 






9009 
ZUUZ 


1 
lU 


12 


7 


DM 1 n Q 1 




Q 


0.43. 10. J 


42. 1 1 .33.3 


23. JO 


-2o0.1 






9009 
ZUUZ 


1 
lU 


1 T 
12 


1 


riN_iU_3_Z 




1 A 


0.43.0y.j 


42. 14.24.2 


24.87 


■2 1/11 

-314.2 






90O9 
ZUUZ 


1 
iU 


12 


1 






1 1 
i i 


0.4Z.4y.D 


42. 13.03. / 


24. 3 


1 C\l 1 
10 /.3 






9009 
zwuz 


1 
iU 


1 9 
12 


-7 








U.44.U0.D 


42.U0.UU.U 


2 J. 02 


-2U2.J 






2002 


10 


13 


Q 







13 


0:44:09.8 


42:06:40.2 


22.99 


-108.0 






2002 


10 


13 


8 


PN.10.3.7 




14 


0:43:31.0 


42:09:25.7 


24.20 


-136.6 






2002 


10 


13 


8 






15 


0:43:39.8 


42:13:46.3 


25.07 


-210.8 






2002 


10 


13 


8 






16 


0:43:37.6 


42:06:58.6 


24.65 


-554.4 


E 


Stream? 


2002 


10 


13 


8 






17 


0:44:59.8 


42:07:44.8 


21.76 


-193.6 






2002 


10 


12 


9 


PN.9.3.4 


1091 


18 


0:45:00.9 


42:08:44.8 


23.53 


-72.2 






2002 


10 


12 


9 






19 


0:44:30.5 


42:08:55.2 


22.36 


-81.9 






2002 


10 


12 


9 


PN.l 0.3.8 


1005 


20 


0:44:51.4 


42:09:34.2 


23.81 


-148.4 


E,R 




2002 


10 


12 


9 







° Probable H II regions (or background galaxies): E — extended in the PN.S data, R — Flux ratio R below cutoff value (see Figure ITTl . other - 
excluded as a FN due to flux . 

^ Sources of non-M3 1 objects: M32 / M32? - Objects associated with and in the vicinity of M32. NGC205 / NGC205? - Objects associated with and 
in the vicinity of NGC205. AndIV - Objects associated with the galaxy Andromeda IV. 2MASXi - Non-M31 source in the 2MASS catalogue. 
MLA93 — Non-M31 source listed in Mevssonnier et al. 1 1993). NS - Objects in the region of the Northern Spur Stream? - Objects referred to in 

iMerrett et alJ l2003l) as being a possible extension o f the Southern Stream . 

= cross-identificatio ns from other catalogues: H 06 - iHallidav et'ai]<20O6h : C89 - ICiardullo etai]<1989ft : HK04 - iHurlev-Keller et aljEooi) (format: 
tablcID); MLA93 - iMevssonnier et alj7l993l) . 



is not much that we can do in terms of detailed dynamical mod- 
eUing. A simple Gaussian fit to the velocity distribution yields a 
mean velocity of —174 ± 6 kms~^ and a velocity dispersion of 
36±5 kms-\The mean velocity is not in particularly good agree- 
ment with the usually-adopted value; the difference seems to arise 
from the statistical coincidence that most of the PNe in M32 reside 
on one side of the galaxy's minor axis (see Figure 1211 . so the mea- 
sured mean velocity contains an element of the system's rotation; 
the somewhat low value for the velocity dispersion is also consis- 
tent with this interpretation. 



7.2 NGC 205 

NGC 205 (M 1 10) is a dwarf elliptical galaxy (dE5) projected ~35' 
from the centre of M31, close to its minor axis. Its major- and 
minor-axis diameters are ^^21. '9 and ~11.'0 re spectively at a po- 
sition angle of 170° Ide Vaucoule urs et all l99l|). It has a mea n ve- 
locity of IINGC205 = -241 ± 3 kms~^ feender et aljl99lh . and 
a mean velocity dispersion of crNGC205 = 46 ± 8 kms~^ dip- 
ping to ^ 20 km s~^ at the very centre of the ga laxy iMateol 19981 
ICarter & Sadler"l990l lsimien & Prugnief200^. 

To identify PNe in NGC 205, we have selected the objects 
within 12.'3 of NGC 205's centre (1.5 times the mean 25.0 mag 
arcsec"'^ B-band radius), and applied velocity cuts at —375 km s~^ 
and —102 kms~^ to seek to isolate the galaxy's members. These 
velocity selections lie at approximately ungc205 ± 3aNGC205, al- 
though the lower bound is placed slightly high to avoid excluding 



43 - 



42 



V 



41 



40 



39 b_l 




2MASXi J0039374+420956 
O 



NGC 205 



Arid VIII 



0.9 0.8 0.7 

RA (hours) 



0.6 



Figure 20. Positions of other galaxies and the the detected emission line ob- 
jects. An object's size and colour reflects its velocity with respect to M3rs 
system velocity: grey objects are receding and black ones are approaching. 
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200 



^ (arcmin) 

Figure 22. Kinematics of emission-line sources in Andromeda IV. The co- 
ordinates are defined relative to the centre of the satellite. Point sizes are 
proportional to velocity with respect to And IV, with grey stars and black 
crosses respectively representing positive and negative values respectively. 



M31 PNe that are seen to have such velocities at slightly larger 
distances from NGC 205. As with M32, the possible overlap at 
these velocities means that there is unavoidably some ambiguity 
in the identification of satellite members, so Table |3| annotates 
those objects that meet the velocity constraints as "NGC205," while 
those that are just spatially coincident are marked as "NGC205?" 
A Gaussian fit to the velocity distribution of likely members yields 
a mean velocity of —235 ± 6 kms~^and a velocity dispersion of 
35 ± 9 kms~^, in agreement with the published values. 

7.3 Andromeda IV 

Andromeda IV is a background dwarf irregular galax y lying a 
projec ted distance of 40' from the centre of M3 1 . Ferguson et alj 
i2OO0h have previously identified a number of emission-line objects 
associated with this galaxy, and we list the cross-identifications in 
the curre nt survey in Tab le|4] The only object that we fail to recover 
from the lperguson et alJ <2000^ study is their "Fl" source, which 
was originally detected in H a, but h as no [O III] cou nterpart in the 
magnitude range of our survey. The lFerguson et aU (2000) objects 
F7 and F8 are located a little away from the galaxy centre, and our 
measurements show that their velocities are consistent with the disk 
of M31, so they do not appear to be members of And IV, and they 
are excluded from the remaining analysis. Table |4| also adds two 
new emission-line sources to the list of objects with positions and 
velocities consistent with membership of And IV. 

The resulting list of seven likely And IV members has a mean 
velocity of 241 ± 10 km s~ ^ , in agreement with the published value 
of 25 6 ± 9 km s ~^, as measured from spectra of only three sources 
iFer guson et all EoOOl) . As Figure l22l illustrates, the velocities do 
not appear to be randomly distributed through the system, but indi- 
cate a line-of-sight component of rotation of ~ 35 kms~^ about a 
centre of RAo ~ 0'M2'"32!3, Deco « 40°34'18'. 

7.4 2MASXi J0039374-I-420956 

Two PN.S emission-line objects that exhibit velocities inconsistent 
with M3 1 , one of which is extended, were found in the vicinity of 
this extended infrared 2MASS source (marked as 2MASXi in Ta- 
ble|3}. The two objects lie at projected radii of 2'.'8 and 24'.'2 from 




-200 - 



(degrees) 



Figure 23. Plot of major-axis coordinate versus line-of-sight velocity, high- 
lighting the postulated location of the And VIII satellite. The small dots 
show the PN.S survey data with Y coordinates between —0.175° and 
—0.450°, the region in which the known And V III candidate PNe lie. The 
open circles are the iHurlev-Keller et all data set within the same 

slice in Y, with the highlighted points indicating the suggested members of 
And VIII. 



the source centre. With only one emission line detected, we can- 
not unequivocally confirm its nature, but if we assume it originates 
from [O III], then this system would be placed at a heliocentric ve- 
locity of ~300 kms^^. As this velocity is similar to that derived 
for And IV, it is possible that both these galaxies belong to the same 
nearby background group. However, further observations would be 
required to confirm that this is the case. 

7.5 MLA93 0953 

On e object with a velocity inc onsistent with M3 1 was also listed in 
the lMevssonnier et alj il993h catalogue of M3I emission-line ob- 
jects. They identified it as either a Wolf-Rayet star or a background 
QSO. Given our velocity measurement, it seems most likely that the 
object is a background QSO with an emission line that happens to 
lie within the [O III] bandpass. 

7.6 Andromeda VIII 

iMorrison et alj J2003h pointed out a coincidence in position and 
velocity of a number of globular clusters and PNe from the 
Hurlev-Keller et al. (2004) survey. These objects appeared to be 
quite well separated in velocity from both the modelled kinemat- 
ics of MSl's thin disk in thi s region and the nearest other PNe, 
leading iMorrison et alj i2003h to suggest that these sources might 
be providing the signature of a new tidally-distorted satellite galaxy 
in M3I, Andromeda VIII. 

However, as Figure l23l illustrates, with our larger sample of 
PNe this gap in velocities between the disk of M31 and And VIII is 
no longer devoid of PNe. The errors on the new velocities are some- 
what larger than those previously published, but not sufficiently so 
to fill in the apparent ~ 50 km s^^ gap in the original smaller data 
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Table 4. Emission-line objects found in the vicinity of And IV, and comparison to the lFerguson et a\\ jlOOOj) data set. 



PN.S 


survey" 










Ferguson et al. (2000)* 


Notes 


ID 


RA 


Dec 


r 


"15007 




ID 






21 10 


0:42:20.3 


40:32:31.8 


174 


25.27 


-386.7 


r 1 




Belovi' flux ratio cutoff Belongs to ^^131 


2111 


0:42:35.3 


40:33:51.0 


44 


24.14 


219.2 






Not found by^e^usoi^^f 


2112 


0:42:21.2 


40:33:50.0 


130 


21.74 


-372.5 


F8 




Belongs to M31 


2113 


0:42:32.2 


40:33:58.8 


20 


22.53 


236.2 


F3 


244 ± 15 


Extended 


2114 


0:42:36.6 


40:34:04.8 


51 


23.47 


213.3 


F2 






2115 


0:42:31.7 


40:34:11.1 


10 


23.13 


253.9 


F4 


250 ± 13 


Extended 


2117 


0:42:32.7 


40:34:17.0 


5 


25.20 


248.8 


F5 






2119 


0:42:30.6 


40:34:47.3 


35 


24.43 


285.7 


F6 


273 ± 19 


Extended 


2128 


0:42:31.9 


40:33:58.7 


20 


24.79 


228.5 






Not found bv Ferguson et al. f2000) 




0:42:37.6 


40:38:12 


241 






Fl 




Not found in PN.S survey 



° ID, RA, Dec, r?i5007 and iihciio frorn Tablelslexcept object Fl, which was not detected, r is the objects radial distance from And IV 

measured in arcseconds^ 

^ ID and velocity froml Ferguson et i200d) Tables 1 and 2. 
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Figure 24. Overlay of the major-axis coordinate versus line-of-sight veloc- 
ity for the PN.S data in the location of And VIII (circles) and the equivalent 
phase-space location on the opposite side of M31 (stars). The PNe previ- 
ously identified as possible And VIII members are highlighted in grey. 



100 




Figure 25. The planetary nebula luminosity function of M31 as derived 
from the current survey. Filled circles are plotted up to the magnitude com- 
pleteness limit; open circles are shown beyond this point. 



tent with membership of M3 1 , it is nearly 0.4 magnitudes brighter 
than any other PN in M3 1 . Unfortunately, it lies outside the region 
mapped by Massev et al. l2002)'s imaging, so we do not have the 
line ratio diagnostic as a further discriminant. However, given the 
very sharp cut-off in the PN luminosity function (see Section[8j, it 
is most unlikely that this object is a PN, so we flag it in the cata- 
logue and exclude it from the following analysis. 



set. Perhaps more tellingly, Figure l24l shows that the equivalent re- 
gion on the opposite side of M3 1 contains a very similar distribu- 
tion of PNe. In fact the equivalent regions in all four quadrants of 
M31 are similarly populated by PNe. We therefore conclude that 
while these PNe may not be members of M3rs thin disk, they do 
probably constitute a normal component of the system's PN distri- 
bution, and it is not necessary to invoke the external explanation of 
a tidally-distorted dwarf companion. 

7.7 Other objects flagged in the PN.S catalogue 

On the basis of the analysis in this section, all objects whose kine- 
matics are inconsistent with M31 have either been identified with 
known satellites, or they have been previously flagged as non-PNe 
on the basis of their extended nature or line ratios (see Section|5}- 
The only other source that we flag as likely not an M3 1 PN is object 
2654. Although this source's kinematics and position are consis- 



8 THE PLANETARY NEBULA LUMINOSITY 
FUNCTION 

Planetary nebulae are observed to follow a well-defined luminosity 
function of the form 

A^(M5007) OC e'^-^°'"''^007 _ g3{Af;oo,-A/5oo7)^ (9) 

JCiarduUo et'ai]|l989l) . with an equivalent formula for apparent 
magnitudes for any individual system. The absolute magnitude of 
the bright-end cutoff, M5007, has a surprisingly universal value 
of —4.5, varying onl y weakly with the host galaxy metallicity 
iCiarduUo e t al.'2002'). Fi gure l25l show s the PNe from this survey, 
with non-members of M3 1 removed, fitted by this function. As the 
figure shows, this functional form fits very well all the way to our 
completeness limit of m5oo7 = 23.75. The apparent magnitude 
of the bright end cutoff, for the sample is msgor = 20.2 ± 0.1. 
This result is in excellent agreement with the value published by 
ICiardullo et ^r\9m of 20.17. 
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Figure 26. The planetary nebula luminosity function at different depro- 
jected disk radii in M3 1 . 



The large number of PNe in our survey means we can divide 
up the planetary nebula luminosity function (PNLF) and examine 
it at different points within the galaxy to test the universality of 
its form. PNLFs for sections at different deprojected disk radii are 
shown in Figure l26l The positional variation in the completeness of 
the survey is clearly shown in this figure: as discussed above, the 
high surface brightness and poor seeing of the central fields mean 
that these data are only complete to msooT ~ 23. For radii between 
0.2° and 1.0°, incompleteness becomes apparent at msooT ~ 24, 
while at the largest radii, for which we used the best seeing condi- 
tions (see Figure|3}, we are close to complete at ms opy ~ 25. 

R ecent obse rvations of PNe in the SMC (.Jacobv & De Marcol 
l2002h and M33 jCiardullo er^l2004 have detected deviations 
from the simple monotonic function of equation j9|. In these ob- 
jects, dips have been observed in the PNLF, whose origins have 
been proposed as resulting from a recent star formation episode. 
The only indication of a dip in the PN.S data set is the point at 
23.2 mag in Figure l25l which lies ~ 2.5a from the model value. 
It is interesting to note from Fi sure l26l that most of the signal re- 
sponsible for this dip comes from the radial range 0?6 < r < 0?8 
where the data dips ~ 2a below the model value at 23.2 mag. This 
is coincident with the star formation ring where the majority of H II 
regions are located discussed in Section|^ suggesting that star for- 
mation may be responsible for this marginally-detected feature as 
well. 

More radically, Ma rigo et alJ J2004) have suggested that the 
brightest PNe all originate from intermediate-mass stars, so that 
the bright end of the PNLF should vary dramatically between 
spheroidal and disk populations. The absence of any significant 
variation in M3rs PNLF with radius, as one goes from the bulge- 
dominated to the disk-dominated parts of the galaxy, implies that 
no such effect is seen in this system. 

As a further test of the universality of the PNLF, and to search 
for the effects of obscuration, we have also divided the data accord- 
ing to whether it comes from the near side (positive Y) or far side 



10 



: 1 ■ ' ' 1 ' ' ' 1 ' ' ' 1 ' ; 
: r < 0.2 , - ! 

* - " ' ^ ^-"""^ - 

■ ** 

/ * 

■ * • CO 

r|M 1 1 1 1 1 1 |T 1 'i 1 ii 


l \ ' ' ' \ ' ' ' \ ' ' ' \ ' I 
': 0.6 < r < 1 .0 : 

i Jm 

■ Q« o o 

■ « ■ 

ri i 1 1 1 1 1 1 1 1 1*1 


■ 0.2 < r < 0.6 

r 1 

: ■ * * : 

- + • ■ 

ip # 

r\ I . . \ , , , I , , ,*\ .- 


: 1.0< r< 2.0 ..i 

*V ■*" \. * * 

- .* J.... » * - 

: Y''* * -. 

^1 *. . t ,- 



20 22 24 26 20 22 24 26 



Figure 27. The planetary nebula luminosity function on the near and far 
side of M31 over different radial ranges. Black stars and the dotted line 
represent the data and fit to the near side of M3 1 (to the north-west), while 
grey circles and the dashed line show the corresponding information for the 
far side (to the south-east). 




Figure 28. An illustration of the location of an Earth based observer with 
respect to the disk of M3 1 . The dust lane is shown as the thick lines in the 
disk plane. On the near side of the disk the majority of light comes from 
behind the dust lane while opposite is true for the far side. Not to scale. 

(negative Y) of the disk. As Figure IZTl illustrates, at most radii the 
PNLFs are statistically indistinguishable. However, once again in 
the region of the ring at 0?6 < r < 1°, we do see a small but sig- 
nificant difference, with the PNLF from the near side systematically 
fainter than that from the far side: values of the best fit for mgoor 
on the near and far side are 20.34 and 20.15 respectively. This 
rather counter-intuitive difference can be qualitatively explained if 
the ring of star formation also contains a modest amount of dust 
obscuration (Walterbos & Kennicutt 1988; de Jong et al. 2005). As 
Figure l28l shows. on the far side of the disk, the line of sight of an 
observer through a highly inclined system like M31 passes through 
more of the galaxy in front of the dust layer than behind, while 
the opposite is true on the near side. Hence on the near side of the 
disk most of the PNe observed are located behind the dust layer, 
dimming the PNLF relative to that seen on the far side. The lack 
of similar differences between other parts of M3 1 suggests that ob- 
scuration is not a major issue across most of the galaxy. 



9 SPATIAL DISTRIBUTION OF PLANETARY NEBULAE 

In order to interpret the kinematics of the PNe in this survey, we 
need to know what stellar population or populations these objects 
are tracing. It is generally assumed that PNe all originate from low- 
to-intermediate mass stars, so they should mainly trace the old stel- 
lar population, but, as we have discussed in Section|8| it has been 
suggested that some fraction of PNe originate from more massive 
stars. As a further test of this possibility, we can compare the dis- 
tribution of PNe with that of the red light from the galaxy, which 
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Figure 29. Surface brightness profile and PN number counts along the ma- 
io r axis of M3 1 . The solid line sh ows the R-band surface brightness profile 
of lWalterbos & Kennicut1<1987|l and the dashed line is an exp onential pro- 
file with a scale length of 0?43 iWalterbos & Kennicutll983) . The nuinber 
density of PNe with msooT < 24 is indicated by closed circles where the 
survey is complete both spatially and to msoor = 24; open circles are used 
elsewhere. 



we would expect to be dominated by relatively low-mass red giants, 
and hence would only be expected to follow the distribution of PNe 
if they originate from similar low-mass systems. To this end, we 
have obtained the radial distribution of the PNe with msoo? < 24 
(the completeness limit of most of the survey) along the major and 
minor axes of M3 1 . 

Along the major axis. Figure l29l shows excellent agreement 
betwe en the number counts of PN e and the R-band photome- 
try of IWalterbos & Kennicutll ^8^). Fitting an exponential func- 
tion to the distribution of PNe beyond the bulge and out to 1.6° 
yields a scale-length of Ra = 0.43 ± 0.02° (5.9 ± 0.4 kpc), 
also i n excellent agreement with the R-band value of 0.43 ± 
0.02° iWalterbos & Kennicutll [l98^ . There are no indications of 
any cut-off in the exponential disk out to the limit of the survey 
at ~ 5Rd. 

Figure l30l shows a similar plot for the minor axis. Once again, 
there is a good match to the published R-band photometry outside 
the central region, including reproducing the bump at F ~ 0.6° 
due to NGC 205. The PNe also follow the published i?^/" profile 
fit to the phot ometr y with a bulge effective radius of i?eff = 0? 10 
(1.4 kpc) jlrwin et alj200l . Even at ~ 10-Rcff > there are no signs 
of an excess over this fit, and hence no evidence for a separate halo 
population even at these very faint surface brightness levels (equiv- 
alent to 26 mag in the R band). 

It thus appears that the PNe trace the same stellar population 
as the old stars seen in the galaxy's red light. They also seem to 
broadly fit a simple two-component model comprising an i?^^^ 
central bulge and an exponential disk out to the very faint sur- 
face brightness levels accessible using such discrete dynamically- 
selected tracers. 
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Figure 30. Surface brightness profile and PN number counts along the mi- 
nor axis of M3 1 . The solid lin e is the R-band surface brightness profiles of 
IWalterbos & Kennicut] ll98i). and the dashed lines is the if.^/"' profile fit 
to the photometrv'bv lh'wh^t^l l2005t) . The number density of PNe with 
"^5007 < 24 is indicated by the filled circles. 
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Figure 31. Variation of 02.5 with distance along the major and mi- 
nor axes. 02 .5 is estimated at positions along the two axes where 
IWalterbos & K ennicutt 1 1987) list a value for the surface brightness, the 
number density of PNe at each location is estimated from the number of 
PNe in a box whose dimensions are equal to the horizontal error bar. 



10 THE LUMINOSITY-SPECIFIC PN DENSITY 

Having established the proportionality between the stellar contin- 
uum emission and the distribution of PNe, we now turn to the con- 
stant of proportionality. This quantity, the luminosity-specific PN 
number density, is usually usually parametrised by the number of 
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PNe in the top 2.5 magnitudes of the PN luminosity function per 
unit bolometric solar luminosity, and is given the symbol 02.5 

The bolometric luminosity of any part of M3 1 can be calcu- 
lated via the standard relations 



Mioi =Mv + BCv 



(10) 



and 

LM31,bol 



= 10 



-OA{Mm 



where BCy is the V-ba nd bolometric correctio n which has a value 
of -0.80 mag for M31 (' Ciardullo et aL 1989). The solar absolute 



is taken to be +4.74 mag 



bolom etric magnitude, MQ,boi, 
I2OOOI) . 

The number density of PNe along the major axis and the V- 
band photometry of Walterbos & Kennicutt ( 1987) [correcting to 
absolute magnitudes using Av = 0.206 ( Schleael et alll99^ and 
a distance modulus of 24.47] give an average value of 02.5 = 
(21 ± 1) X (lO'' l/0,i,o!)~^- FigurelsTlshows, outside the very 
uncertain central bin there is no strong evidence for any gradient 
in the value of 02.5 with position, as we would expect on the basis 
of the close proportionality determined in the previous section. The 
same calculation on the minor axis, where the light profile is domi- 
nated by the bulge, gives 02.5 = (15 ± 2) x (10^ LQ^boi)~^, once 
again with no evidence for a gradient (see Figure IsTl . This value 
agrees very well with the calculation by Ciardullo et al. {J3W^ 
for the PNe in MSl's bulge, from which they derived a value of 
02.5 = 16.3 X (lO'J L0,boi)~'- 

The modest difference between the mean values for the minor 
and major axes is consistent with the existing evidence that Q2.5 is 
correlated with galaxy colour, with red ellipticals being poorer in 
PNe per unit lu minosity than blue spiral galaxies <PeimberJl990l: 
lHuietalll99l : here for the first time we are seeing evidence for 
such a difference between spheroidal and disk components within 
a single galaxy. Fortunately, such a mild deficit in the number of 
PNe from the bulge component relative to those from the disk will 
have little impact on the dynamical modelling of these data, since 
it merely changes the effective normalisation of the number den- 
sity of tracers in these two discrete components. We can therefore 
proceed to make a preliminary study of the PN kinematics in M3 1 . 



11 PLANETARY NEBULA KINEMATICS 

In Figure l32l we show a plot of PN velocity against major axis po- 
sition. A number of fairly subtle dynamical features are apparent in 
this plot. For example, many of the satellite galaxies identified in 
Section|7|(shown here as circles) stand out quite clearly. Similarly, 
the PNe identified by iMerrett et alJ i2003t) as a possible continu- 
ation of the merging Southern Stream of stars hbata et al]l200lL 
|200# .' McConnachie et al. "2003"), shown as triangles, are apparent. 
Interestin gly, th e other strong photometric structure, the Northern 
Spur ^Ferguson et al . 2002), does not stand out in the plot. The 
kinematics of PNe in this region (shown as squares in Figure l32l 
are indistinguishable from those of the disk, despite the Northern 
Spur's location some distance from the plane. It would thus appear 
that this feature is related to the disk, perhaps indicative of a warp. 
After these minor features have been excluded. Figure l32l shows 
that the kinematics are dominated by rotation at large radii, but that 
at small radii there is a larger random component, attributable to a 
bulge population. This result indicates that the majority of the PNe 
form a relatively simple system of rotationally-supported disk and 
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Figure 32. Scatter diagram of position along the major axis of M3 1 versus 
velocity for all the survey PNe. The highlighted objects ai'e those in the 
region of satellite galaxies (circles; see Sectional; those in the region of the 
Northern Spur (squares); and those identified as forming a continuation of 
the Southern Stream (triangles). 



randomly-supported bulge, in agreement with the conclusion based 
on photometry in Section|9| 

A full dynamical model fit to the entire dataset is beyond the 
scope of this paper, but we can obtain some interesting insights 
by considering the PNe that lie close to MSl's major axis. For a 
highly-inclined galaxy like M3 1 , the observable kinematics in this 
region are dominated by the rotational motion, so we can essen- 
tially isolate this component of velocity. Accordingly, we have se- 
lected a sample of 690 PNe with |y| < 0?04, close to the major 
axis. Figure l33l shows a plot of this subsample's line-of-sight veloc- 
ities versus major axis position, together with the mean velocities 
derived by radially binning these PNe into groups of 50 objects. 
A system velocity of —309 km s^^ has been adopted as discussed 
previously. 

The distribution of PN velocities in Figure |33]has an upper en- 
velope that matches the rotation curves derived from various tracers 
(also shown in this figure) reasonably closely at large radii. At small 
radii, the PNe show a wide spread in velocities, indicative of a con- 
tribution from the bulge. It is also interesting to note the dip in the 
mean rotation speed of PNe at ~ 0?3. This feature, which shows 
up even more clearly in the CO rotation curve, has been attributed 
bv lLoinard et alj j 19951) to non-circular motions associated with a 
barred bulge. Its appearance in the PNe kinematics confirms that 
they are dominated by disk objects even at this fairly small radius, 
but also serves as a warning that large non-circular motions may 
complicate the observed dynamics at these radii. At larger radii, the 
mean velocities of the PNe trace the rotation curve relatively well, 
except that the mean velocity lies systematically below the circular 
speed. We will return to the origins of this asymmetric drift in the 
mean velocities of PNe in Section fl 1.21 but first we look in more 
detail at their random motions. 
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Figure 33. Plot of position versus velocity for PNe close to the major axis of 
M31 < 0?04). The data from both sides of the galaxy have been com- 
bined, and the velocities have been corrected for the galaxy's inclination. 
The grey-scale shows the density of PNe in this phase-space projection. 
The circles show the mean velocities of these data in a number of bins - 
the horizontal error bars show the extent of each bin, and the vertical bars 
the un certainty in m ean velocity. The dotted line shows M3 1 's H I rotation 
curve iBraunll99 l'): the dashed line shows the rotation curve derived from 
CO measurements iLoinard et al. 1995); and the solid line shows a rotation 
curve derived from H II region velocities lKenil989l) . 
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Figure 34. Velocity dispersion profile for the major-axis PN sample shown 
in Figure l33l Points with error bars show the binned PN dispersion profile, 
corrected in quadrature for the uncertainty in the PNe velocities. The solid 
grey line is the inner stellar velocity dispersion profile of McElrov 1 1983). 
The curves are for a number of models fit to the PN data with a minimi- 
sation: the grey dot-dash line shows the empirical fit of equation 1121 : the 
solid black line shows a simple exponential de crease; the dotted black U ne is 
for a flaring disk model with parameters from lBrinks & BurtonHl984l) ; and 
the dashed black line shows a flaring disk model with parameters matched 
to the PN velocity dispersion profile. 



11.1 Velocity dispersion 

Figure |34] shows the observed velocity dispersion as a function of 
position along the major axis obtained from the PNe in Figure l^ 
binned radially into groups of 50 objects. Near the centre, the cal- 
culated PN vel ocity dispersions agree well with the declining pro- 
file derived bv lMcElrovl il983h from stellar absorption lines (also 
shown in the figure), providing further confidence in the PN results. 
However, the stellar absorption-line data only reach to ~ 0? 2, while 
the PN profile goes out to ~ 2° (27.4 kpc). It is in these previously- 
unexplored outer regions that the velocity dispersion starts to be- 
have a little strangely. First, there is a dip in the observed disper- 
sions at ~ 0.3°, suggesting that the bar's influence shows up in the 
dispersion profile as well as the mean rotational velocity. Second, 
the dispersion profile seems to flatten out at larger radii to a con- 
stant value of ~ 55 kms~^. Ignoring the bar-related dip, a simple 
empirical fit to the large-scale profile, plotted as a grey dot-dash 
line in Figure l34l is given by 



cr^ = aA + CTB exp 



(--) 

\ RdJ 



(12) 



with a A = 47 kms"\ as = 84 kms"^ and Ra = 0?43. Other 
recent stellar kinematic data from large radii also reveal significant 
dispersions: Ibata et al- (2005) find velocity dispersions of RGB 
stars of ~ 30 kms~^ from small fields at even larger distances 
along the major axis. 

These results are not consistent with the expectations of the 
simplest disk models. By applying the tensor virial theorem to a 
one-dimensional sheet approximation to the vertical distribution of 



Figure 35. Schematic diagram of the flaring disk of M31. The disk scale 
height is constant out to a radius it then increases linearly with radius 
to -Rout . at which point it becomes constant again. 



stars in a disk, one can show that the vertical velocity dispersion, 
a^, should be given by an equation of the form 

2 



o"j. = 27rGE2o 



(13) 



where zp is the scale heigh t and E the surface density of the disk 
jBinnev & Merrifieldll99a. p.727). If the surface density of stars 
drops off exponentially with radius with a scale-length Rd, and the 
scale-height stays constant with radius, then should drop expo- 
nentially to zero with a scale-length of 2Rd. If the velocity ellipsoid 
remains roughly constant in shape, then the other components, in- 
cluding the observed <j) component, should decrease in the same 
way. As the solid line in Figure l34l shows. this is not what we find 
in M31. 

So, where have we gone wrong? Well, it is known that the 
outer parts of MSl's disk are not a simple constant-thickness flat 
sheet . The outer H I disk is s ignific antly warped ( Brinks & Shan^ 
'1984*), and 'Brin ks & BurtonI ^984*) have suggested that the thick- 
ness of the H I layer also flares to greater scale-heights at large radii. 
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Their model for the flaring, shown schematically in Figure |35| as- 
sumes that the Hi scale-height is constant at Zin = 0?01 (136 pc) 
out to 1?33 (18 kpc), then increases linearly out to _Rout = 2? 54 
(34.5 kpc) where Zout = 0? 14 (1.9 kpc). Be yond this radius the H I 
becomes undetectable, but llbata et aljilOO S') have assumed the disk 
scale height becomes constant from this point onwards to account 
for their stellar observations. 

If the scale-height increases in this way, then equation <13> 
implies that, for a given mass density, the vertical velocity dis- 
persion also increases. If we insert the above prescription for the 
scale-height into this equation, but continue to assume that the disk 
density drops exponentially with radius, then we find the velocity 
dispersion profile shown as a dotted line in Figure l34l Although 
a step in the right direction, the fit to the observed profile is still 
not very good. Of course, there is no reason why the stellar disk, 
having been built up over an extended period and undergone sec- 
ular heating since it formed, should flare in the same way as the 
H I distribution. Indeed, it is often assumed that the stellar scale- 
heights of disk galaxies remain constant with radius. However, re- 
cent analysis has shown that existing photometric data do not rule 
out si gnificant variation in scale-height with radius ( Naravan & Jog 
l2002h . Further, if a disk comprises two components of different 
constant scale-heights, with the "thick disk" having a longer scale- 
length than the "thin disk," then one would expect a variation with 
radius in their combined effective scale-height similar to this sim- 
ple flaring model. Whether such components should be combined 
or should be considered as discrete entities is largely a matter of 
semantics unless one can establish some other feature that distin- 
guishes their constituent stars. If we use this toy varying-scale- 
height model, but move the point at which the flaring starts to 
Rin = 0? 85 (1 1.6 kpc) while leaving all the other parameters unal- 
tered, we obtain the dashed line in Figure l34l in respectable agree- 
ment with the data. Clearly, by altering the way in which the layer 
flares, or by varying the mass-to-light ratio of the disk with radius, 
we could improve the fit still further, but such detailed modelling 
goes beyond what is justified by the simplifying approximations of 
a one-dimensional sheet of stars in hydrostatic equilibrium. 

11.1.1 Variation with magnitude 

This analysis of the velocity dispersion allows us to make a fur- 
ther search for any systematic variation in the properties of PNe 
with luminosity. Such a test is motivated by knowledge of the Solar 
neighbourhood of the Milky Way, where we know that the tangen- 
tial components of the velocity dispersions of different types of star 
vary by more than a factor of two, with the youngest, most massive 
stars having the lowest dispersions (Binnev & Merrifield 1998). If 
the PNe in M3 1 that populate different parts of the PNLF had pro- 
genitor stars with very different masses, then we might expect sim- 
ilar variations in their velocity dispe rsions. Not only wo uld such 
a variation fit with the suggestion of 'Marieo et al." ("2004") that the 
bright end of the PNLF is populated by mo re massive stars, but i t 
would also resonate with the recent claim bv lSambhus et alJ i2006h 
that the PNe in NGC 4697 can be divided on the basis of their lu- 
minosities into two discrete populations with distinct kinematics. 

To test for such an effect, we have extracted from the major- 
axis subsample those objects for which X > 0.3°: since the global 
velocity dispersion is approximately constant at these radii (see 
Figure l34l . we can combine data across this range to form a mean- 
ingful average. Figure 1361 shows the velocity dispersions deter- 
mined from these data combined into half-magnitude bins. Clearly, 
there is no evidence for an decrease in dispersion for brightest ob- 
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Figure 36. Variation in velocity dispersion with apparent magnitude for the 
major-axis data from the disk region at Jf > 0.3°. 



jects, adding further weight to the conclusion that PNe of all mag- 
nitudes form a relatively homogeneous population of old stars. 



11.2 Asymmetric drift 

The presence of significant random velocities at all radii means 
that the support against gravitational collapse does not come en- 
tirely from circular motions even at large distances from the centre 
of M3 1 . We would therefore expect the observed mean rotational 
streaming of the stars, v^, to lag behind the circular speed, Vc. As 
noted above, such an asymmetric drift is seen, but having measured 
the velocity dispersion profile we can now quantify whether the 
predicted velocity difference is dynamically consistent with the ob- 
servations. 

The relationship between and Vc at an y radius R is given by 
the ap propriate axisymmetric Jeans equation jBinnev & Tremaina 
119871 equation 4-33), which we can rewrite as 



V4, = 



„2 J ^ Rl _|_ ^ 



(14) 



''^^"R^ u OR ■ " dz 

where we have adopted a cylindrical coordinate system, {R, 0, z}, 
and measure stellar velocities {vR,v,iy,Vz} and number density 
iy{R, z) in this reference frame. To solve this equation, we have to 
make some further assumptions. First, the final term, which arises 
from any tilt in the velocity ellipsoid, is unobservable. However, 
for all plausible dynamical models it remains small and has lit- 
tle impact on an analysis of this kind (Gerssen et al. 1997), so we 
can safely neglect it. Second, we must assume something about the 
shape of the velocity ellipsoid of random motions. The simplest as- 
sumption is that it remains constant with position, so that we can 
write 



(15) 



where a is a constant. In the solar neighbourhood o f the Milky 
Way, a is found to be ~ 2 iBinnev & MerrifielJl99^ . and in the 
absence of any other constraint we set it to this value at all radii in 
M31 (varying a by ±0.5 makes little difference to the final result). 
If we make this approximation, use the H I data to obtain a simple 
linear approximation to the rotation curve Vc{R), take the number 
counts of PNe in Section|9|to give the exponential disk distribution 
for ^{R), and use the smooth approximation to (j^{R) of equa- 
tion J12> . then we can substitute these quantities into equation <14t 
to predict the value of tJ^. 

To compare this value to the observed mean streaming of stars, 
we must make a couple of corrections to the model. The first is 
due to the inclination of the galaxy, but M31 is sufficiently close 
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Figure 37. Predicted and observed rotational motions of PNe in M3 1 . The 
filled circles show the observed streaming of stars, as in Figure l33l The 
crosses give MSl's H I rotation curve as derived byjBraunji 1991), with the 
straight line fit to the data used in this analysis shown by the dotted line. 
The solid line shows the predicted rotational velocity obtained by solving 
equation 1141 . 



to edge-on that this effect is relatively minor and has already been 
accounted for in the data of Figure l33l The second effect is more 
significant: although we have selected a sample of PNe whose po- 
sitions lie close to the major axis of M31, the cut at |F| < 0?04 
does sample a range of azimuths, particularly at small radii. We 
have allowed for this effect by averaging the line-of-sight compo- 
nent ofv^ in the model over the range of azimuths admitted at each 
position X along the major axis. 

The net result of this calculation is shown in Figure |37l At 
small radii, the detailed match between model and observation is 
not perfect, but discrepancies here are not surprising since many 
of the PNe in this region will belong to the bulge rather than the 
disk as modelled here; even PNe from the disk will have their 
kinematics affected by the bar, which also has not been accounted 
for in this simple axisymmetric picture. At larger radii, however, 
the agreement is remarkably good. It should be stressed that this 
model has not been fitted to the data: the amplitude of the asym- 
metric drift between and Vc is entirely fixed by the Jeans equa- 
tion once we have determined the other quantities in this equation 
as described. Note also that although we have constructed this dy- 
namical model by assuming the scale-height of MSl's disk pop- 
ulation increases with radius, the absence of in equation <14t 
means that the agreement is not dependent on this specific model. 
All that is required is that the radial and tangential kinematics are 
coupled via equation <15t : if is not coupled to these motions, 
then any scale-height can be accommodated. The generic agree- 
ment between model and data therefore provides a healthy consis- 
tency check on this analysis, and leaves us with a relatively simple 
self-consistent picture of the dynamics of MSTs surprisingly warm 
disk. 



12 CONCLUSIONS 

Using the Planetary Nebula Spectrograph, we have surveyed a wide 
area over the Andromeda Galaxy. Some 3300 emission-line objects 
have been detected and their velocities measured. Cross-checks 
against other smaller data sets have confirmed the accuracy of these 
measurements, and we have also analysed the detected objects in 
order to separate the PNe from other emission-line objects. Al- 
though this paper is primarily intended just to present the data, we 
can draw a number of intriguing conclusions from even the initial 
analysis of these observations. Specifically: 

• A number of M3rs satellites are detected in the survey, but 
the existence of Andromeda VIII is brought into question by this 
larger data set. 

• The luminosity function is well reproduced by the usual sim- 
ple functional form, with no strong evidence for spatial variations 
that cannot be attributed to modest amounts of obscuration. It has 
been suggested that the bright end of the luminosity function is 
dominated by relatively massive stars, so one might expect the 
PNLF to vary significantly between the spheroid and disk com- 
ponents. The absence of such effects suggests that all PNe in M3I 
arise from old low-mass stars. 

• The number counts of PNe matches that of the R-band pho- 
tometry, further indicating that both arise from the same old popu- 
lation. With these dynamically-selected number counts of discrete 
objects, we can also go much deeper than the conventional pho- 
tometry of M31. There is no evidence for a cut-off in its exponen- 
tial disk population to beyond four scale-lengths, and there are no 
signs of a halo population that is distinct from the bulge out to ten 
effective bulge radii. 

• Somewhat unexpectedly, the velocity dispersion of M3 1 does 
not drop to zero at large radii, and there is a significant asymmetric 
drift between the rotational motion and the local circular speed at 
all radii. These data can be self-consistently reproduced by a simple 
disk model which is partially supported by random motions even at 
large radii. If the shape of the velocity ellipsoid remains approxi- 
mately constant, then M3rs disk must increase in scale-height at 
large radii or its mass-to-light ratio must increase significantly. 

• There is no evidence for any systematic variation in disk PN 
velocity dispersion with magnitude, again suggesting that PNe of 
all luminosities arise from a relatively homogeneous old stellar 
population. 

Clearly, this quick pass through the database has only 
scratched the surface of what can be done with a kinematic sur- 
vey of this size. We look forward to more detailed dynamical mod- 
elling of the full three-dimensional {X, Y, v} phase-space data set, 
as well as tagging the PNe with extra dimensions such as their mag- 
nitudes and line ratios. Such analysis will allow us to dig deeper 
into this rich parameter space in order to further understand the dy- 
namical properties of our nearest large neighbour and its planetary 
nebula population. 
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